Ground-state properties of rutile: electron-correlation effects 
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Electron-correlation effects on cohesive energy, lattice constant and bulk compressibility of rutile 
are calculated using an ab-initio scheme. A competition between the two groups of partially covalent 
J_h ' Ti-O bonds is the reason that the correlation energy does not change linearly with deviations from the 
Oh. equilibrium geometry, but is dominated by quadratic terms instead. As a consequence, the Hartree- 
■ Fock lattice constants are close to the experimental ones, while the compressibility is strongly 
. ' renormalized by electronic correlations. 
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i 1! I. INTRODUCTION 

• i— I 

<-> : 

Although transition-metal oxides are one of the most 
' interesting classes of solids, relatively little theoretical 
-y . work aiming at a microscopic understanding of electron- 
Q ' correlation effects in these systems is available so far. 

, This is not surprising: high numerical effort is already re- 
£2 ' quired here for accurately describing ground-state prop- 
erties at the Hartree-Fock (HFYirjndepcndent-particle 
' ■ level. Recently, results on NiOoH and on the rutile 
TiC>2 crystaocl were published. Correlation effects in 
g , Ti02 have beeustudied only implicitly, at the density- 
q ' functional leveloQ. 

I— ~~ ', Rutile is one of the experimentally found modifica- 
| ■ tions of Ti02- The investigation of the relative stability 
^ | compared to other phases such as anaiase, brookite or 
• Ti02(B) is another interesting subject!]. Rutile is well 
' studied by experimentalists; the structure was precisely 
t-H , determined by using X-ray and neutron diffraction (for 
a discussion of the experimental results see Ref. [| and 
references therein). The experimental lattice constant is 
close to values obtained in ab-initio self-consistent field 
(SCF) calculations^, whereas the difference between SCF 
results and experiment for the compressibility is large. 
This contradiction makes the task of studying electron 
correlations in the rutile crystal interesting. 
*^ ' From the methodological point of view the rutile crys- 
Ch , tal is a next logical step in the application of a corre- 
lation treatment within the so-called "scheme of local 
increments"™]. The validity of this approach was suc- 
cessfully tested for covalently bonded solids like diamond, 
graphite and many typical scmiconductorsBiii For 
purely ionic crystals the scheme works also welL|Csee,|Eef- 
erences for MgO, CaO, NiO and alkali halidesyO'BB). 
II- VI semiconductors have been investigated as, an ex- 
ample of partly ionic, partly covalent crystalscZI. With 
respect to this, rutile is another excellent object for test- 
ing. 

In the present short contribution we report on ab-initio 
correlation calculations for the cohesive energy, lattice 
constant and compressibility of rutile. The paper is or- 
ganized as follows. In the second chapter we outline the 
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computational method and describe our results. We also 
discuss individual correlation-energy increments and ad- 
dress the question of their transferability when approxi- 
mating the infinite crystal by embedded clusters. A dis- 
cussion and a short summary are presented in the last 
chapter. 

II. COMPUTATIONAL METHOD 

The rutile structure is tetragonal (nonsymmorphic 
space group PA^/mnm) with two titanium atoms and 
four oxygens per primitive unit cell (see Fig. |l|) . Ti atoms 
are located at positions (0,0,0) and (1/2,1/2,1/2) and 
oxygens at positions (±2;, ±x, 0) and (1/2 ±x, 1/2 =Fx, Q-L 
The lattice constants extrapolated to zero temperatureE£l 
are a = 4.592 A and c = 2.958 A, the dimensionless coor- 
dinate x = 0.3048. The experimental value for 9 is 98.8°. 
The bonding and charge distribution in the rutile crystal 
have been discussed in Ref. [|. The SCF Mulliken popu- 
lation analysisau shows that the excess negative charge 
on oxygen is about — 1.4e. On Ti, in addition to the 
closed (but easily deformable) argon-like core the popu- 
lation analysis gives 1.2 valence electrons in the d shell, 
whereas the 4s shell is not occupiecocl. 

After this preliminary discussion, we briefly review our 
theoretical approach for including electron correlation on 
top of a crystal SCF calculation. The main idea is to ex- 
pand the total correlation energy, E corri of the crystal 
in terms of local correlation-energy incre ment s. Details 
have been described in previous papersailiUla. Essen- 
tially, we use a Bcthc-Goldstone-like hierarchy of the type 

E corr =^2e(A) + ^J2 Ae ( AB )+ (!) 

A A,B 

1 £ Ae(ABC) + ... 

A,B,C 

where A, B, C, . . . denote groups of occupied localized or- 
bitals on atoms A, B, C, . . ., respectively (we use the same 
symbol for an atom and for the group of the localized 
orbitals on A). The quantity s(A) (a one-body incre- 
ment) denotes the (local) correlation energy of the crys- 
tal in which only the orbitals of group A are correlated. 
The two-body increment Ae(AB) is defined as the non- 
additive part of correlations arising from simultaneously 
correlating the groups A and B: 

Ae(AB) = e(AB) - e(A) ~ e(B) (2) 

Similarly three-body increments can be introduced as fol- 
lows: 

Ae{ABC) = e(ABC) - Ae(AB) - (3) 
Ae(AC) - Ae(BC) - e(A) - e{B) - e{C) 

Even if one should want to treat these formulae as 
purely phenomcnological ones, still the expansion (if car- 
ried out to infinity) is exact which is clear upon in- 
spection. (For a more sophisticated derivation, cf. Ref. 
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|l9| .) For practical applications, however, some approxi- 
mations are inevitable. The first approximation consists 
in truncating .the. infinite series. The data for various 
systemsErHEZln'o1l3 show that the expansion is quickly 
convergent. The three-body increments were found to 
be almost negligible (for the bulk of practical applica- 
tions). Moreover, two-body increments were found to 
decay rapidly for larger distances. The second (more se- 
rious) approximation for the incremental expansion con- 
sists of replacing the infinite crystal with finite embed- 
ded clusters for the purpose of determining individual 
increments. This approximation uses the local nature of 
electron correlationE2l and was found to work well. Note 
that we restrict the local cluster treatment to correlation 
effects only — SCF interactions are long-range, and a cal- 
culation involving the whole (infinite) lattice is manda- 
tory heEe,,-This is possible with the CRYSTAL program 
packageOM 

Applying the combined scheme just described to TiC>2, 
one obtains the total energy functional 

E tot (a,c,x) = E sc f(a,c,x) + E corr (a, c, x) (4) 

as a function of the lattice parameters a, c,iorri, T2, 6 
(in internal coordinates). The crystal compressibility B, 
at zero temperature, is defined as 

d 2 E (V) 

= V —^- (5) 

with V denoting the volume of the unit cell and Eq(V) 
the conditional minimum of _E tot (a, c, x) for fixed constant 
volume. Using V — a 2 c and expressing E tot (a, V/a 2 , x) 
as a function of two independent variables a, x, the com- 
pressibility B can be obtained. This method is usually 
not applied to noncubic crystals, but instead the mini- 
mum position at do, eg, xq is approached by a conjugated 
gradient technique. This reduces the number of data 
points, and the compressibility can only be obtained then 
by applying an empirical ansatz such as the Murnaghan 
equation of state (cf. e.g. Ref. |j). Another possibility 
is to define an artificial isotropic compressibilityu Bjso 
where one (incorrectly) assumes that under compression 
x remains constant and a and c scale isotropically. 



A. Basis sets 

SCF calculations: For the titanium atom, a relativislic 
energy-consistent 12- valence-electron pseudopotentialu, 
together with a 411/411/41 basis seta was used. Fps oxy- 
gen we chose the basis set given by Causa et a£j. To 
calculate the energy of the free atoms, diffuse functions 
cannot be omitted (let us remind J_b at they must be omit- 
ted in CRYSTAL calculationsBB) . Thus we awnropri- 
ately supplemented the basis sets for this purposeEa. Our 
SCF cohesive energy is smaller compared to that of Ref. 
probably just because of the lack of diffuse functions 
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for the free Ti atom in Ref. [4|. Of course, this omission 
affected only the cohesive energy. 

Correlation calculations: For oxygen we used Dun- 
ning's correlation-consistent augmented valence triple 
zeta [5s4p3<i2/] basis setn3. For Ti, our starting point was 
a [6s5p3d] basis set optimized for the pseudopotential of 
Ref. ^3]. We decontracted one p and one d function, end- 
ing up with a [6s6p4<i] set which was supplemented with 
2 fig polarization exponents optimized in CCSD calcu- 
lations for the free ground-state atom (/ exponents are 
2.45, 0.766, and the g-exponent is 2.132). 

B. SCF calculations 

With the CRYSTAL95 program packageEE, SCF en- 
ergies for 100 different sets of randomly chosen a, c, x 
values were calculated. A region around the minimum 
was covered more densely. The computational results 
were fitted by a 3- variable polynomial of third order plus 
a few quartic terms which turned out to be significant. 
(Typically about twenty constants to fit.) 

Some final results of SCF calculations are collected in 
Table I. We recover at the SCF level about 57 % (834 
mH) of the experimental cohesive energy. The lattice 
constant is already close to experiment, but the bulk 
modulus turns out to be substantially too high. These 
results indicate the importance of taking into account 
electron correlation. 



C. Correlation calculations 

For the correlation calculations, we applied the 
coupled-cluster approach with single and double exci- 
tations (CGSD)EZI and included perturbatively triples 
(CCSD(T)jcp "° implemented in the program package 
MOLPROE3E9. Within the incremental scheme, we stud- 
ied one-, two- and three-atom clusters. (For the three- 
atom clusters, we performed CCSD calculations only.) 
The clusters were embedded in a large slab of Madelung 
point charges (7x7x9 unit cells, charges pk4_aiid -2, 
respectively), similarly as described in Refs.EfEjO'EEl. An 
exception are the Ti ions nearest to. the cluster atoms: 
in this case, Ti 4+ pseudopotentialsEl were used instead 
of the bare point charges, in order to simulate the Pauli 
repulsion on the O 2- electrons of the inner cluster. The 
not fully ionic character of the system turns out to be 
the basic difficulty. When doing computations on finite 
clusters, only an integer electron number is allowed and 
the atomic populations obtained from the CRYSTAL cal- 
culation cannot easily be reproduced as is the case for 
perfectly ionic systems. We are forced to assign to the 
clusters considered the same total electron charges as in 
hypothetical purely ionic Ti 4+ , 2 ~ rutile. Of course, 
this is no issue in large clusters, but may be critical for 
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small ones. In order to control the quality of this ap- 
proximation, we compare the results for the individual 
increments taken from clusters with one, two, and three 
explicitly described ions (see Table II). 

The bulk of our correlation-energy calculations (whose 
results appear in Table I) was done for embedded clusters 
with one and two explicitly treated atoms. The maxi- 
mum distance between the atoms within these clusters 
was up to 8 atomic units. These clusters (altogether 13 
in number) provide us with the most significant one-body 
and two-body contributions to E corr . For 5 different ge- 
ometries, all these increments were calculated. For 15 
other geometries, the 8 most important increments were 
explicitly calculated whereas the remaining 5 least im- 
portant (contributions below 1 mH ) were obtained by 
interpolation. To estimate the importance of three-body 
contributions, we also calculated the biggest three-body 
O — Ti — O increments for one geometry, and they were 
found to be of the order of 1 mH. 

A very important question is the transferability of 
these increments from one type of cluster to the other (see 
Table II). In every case, our increments refer to localized 
orbital groups which can be formally attributed to O 2- 
and Ti 4+ ions, respectively, but due to partially cova- 
lent bonding in TiC>2, especially the former orbital group 
becomes more delocalized when going to larger clusters. 
Specifically, we see that the oxygen correlation energy 
increases in magnitude when taken from a cluster with 
three explicitly described ions (one titanium ion, two oxy- 
gen ions, to be consistent with the formula unit) instead 
of out of a cluster where only the oxygen ion has basis 
functions. We explain this with a better description of 
the diffuse tail of the oxygen charge, as the charge is al- 
lowed to flow to the titanium ion. The more diffuse ion 
has lower-lying excitations which explains the change of 
the correlation energy. In the same way, we can compare 
the titanium correlation energy of a cluster where only 
one titanium ion has basis functions to that of a cluster 
with three explicitly described ions. We find that the ti- 
tanium correlation energy is reduced in magnitude. This 
is due to the oxygen charge which has flown to the tita- 
nium ion and has led to a d occupancy. Excitations of the 
3s and 3p electrons into d orbitals are now (partially) for- 
bidden, i.e., we have an exclusion effect. In total, both ef- 
fects almost cancel for the one-body increments (|Ae(0) 
increases by 5 mH, to be multiplied with a weight fac- 
tor of 4, (\Ae(Ti)\ decreases by 11 mH, to be multiplied 
with 2). For the two-body increments, a partial cancel- 
lation with increasing cluster size takes place, too: 0-0 
increments are enhanced by ~4 mH when bridging Ti 
atoms are taken into account, while Ti-0 increments are 
reduced in magnitude by ^3 mH in the Ti02 unito. 

As a whole, we find that the transferability is reason- 
able but rather poor in comparison to purely ionic crys- 
tals where cluster charges are strictly confined (see Table 
II). This is not unexpected as the perfect-ion approxima- 
tion is broken in a different way for each of the clusters 
studied. To overcome this transferability problem, the 
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only way would be to further improve the cluster sur- 
roundings by using more than three explicitly described 
ions with high quality basis set (ideally in multiples of the 
formula unit), which is presently not possible because of 
the steeply increasing computational effort. 

The correlation energies obtained this way were fitted 
for 20 geometries by a second-order polynomial (ten pa- 
rameters), analogously as previously described for SCF 
energies. Adding the two resulting formulae, we obtained 
Etot (a>, c, x) so that we could perform the calculation of 
physical constants which include correlation effects (see 
Table 1). 

III. DISCUSSION AND SUMMARY 

As SCF and SCF+CCSD(T) lattice constants are not 
very different and both rather close to the experimental 
values we conclude that electron correlations do not influ- 
ence the crystal geometry in a major way. The technical 
explanation is provided by studying the coefficients of 
the analytic (fitted) formula of E corr (ri, r 2 , #), in atomic 
units, at the CCSD(T) level: 

E CO rr(ri,r 2 ,9) « -2.1241 - 0.0121Ari (6) 
-0.0067Ar 2 - 0.0322A6* 
+0.2570(Ari) 2 - 1.5730(Ar 2 ) 2 - 1.2490(A6») 2 
-3.0459AnAr 2 + 1.358OAnA0 - 0.0875Ar 2 A# 

The normalized dimensionlcss expansion variables are de- 
fined as Ari = (ri — ricOAiexp where rio is the SCF 
equilibrium value of the internal variable r\ and r\ exp 
is the corresponding experimental value. Analogous ex- 
pressions hold for Ar 2 and for A8. The coefficients of 
the linear terms are much smaller than the coefficients of 
the quadratic terms. This explains why the compressibil- 
ity is strongly renormalized by correlations as it depends 
primarily on the coefficients of the quadratic term. 

To understand this situation it is necessary to study 
how individual increments change when the lattice con- 
stants change. To gain a qualitative understanding it is 
enough to study the biggest contributions: the one-body 
increment for oxygen and the two-body nn Ti-0 incre- 
ments. (Note that the one-body increment for Ti can be 
considered as constant, to a good approximation). In the 
following we will describe several competing mechanisms 
which cause linear changes of the bulk correlation energy 
to be small near the SCF minimum. 

Let us start with oxygen (one-body increment for oxy- 
gen) . Its change as a function of the lattice constant is in 
part of intraatomic origin and in part of electrostatic ori- 
gin. First let us consider the intra atomic part. As found 
and explained in earlier workoEjo, excitations cost less 
energy when the lattice constant (and the volume of the 
quantum well enclosing the oxygen ion) increases; in that 
case, the magnitude of the correlation energy increases, 
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too. There is an opposite trend connected with the per- 
manent electric field E at the oxygen site. As the lat- 
tice expands, the Madelung field decreases, and the in- 
fluence of correlation on the static polarization of the oxy- 
gen ions, \{otsCF — oi C orr)E 2 (with polarizabilities asc'F 
and a corr at SCF and correlated levels, respectively), de- 
creases in parallel. (Note that due to symmetry there is 
no static electric field at the site of the Ti ions). A compe- 
tition between the two effects just described apparently 
leads to a relatively small change in the total oxygen cor- 
relation energy as a function of the lattice constant. 

Thus, the most important source of variability are two- 
body increments. For the two-body nn Ti-0 increments, 
we find two competing effects again. To fix the attention 
let us provide two greatly simplified formulae (atomic 
units): 

Ae api « -0.027 + 0.083Ar! - 0.039Ar 2 - 0.036A6> • • • (7) 
Ae equa « -0.026 - 0.046Ar! + 0.018Ar 2 + 0.040A6> • • • 

where Ae api and Ae equa are two-body nn Ti-0 incre- 
ments for apical and equatorial oxygen ions, respectively. 
The large quadratic terms in the above formulae were 
not shown as they are irrelevant to the following argu- 
ments. Let us remember that the multiplying weight 
factors (per unit cell) for Ae ap i and Ae equa are 4 and 
8, respectively. Due to the opposite signs and due to the 
fact that the different linear term amplitudes are roughly 
in 1:2 proportion we may conclude that a small lattice 
distortion can change the individual Ae but not the sum 
4Ae a pi + 8As equa which is roughly constant (but only 
in linear approximation). Thus, 4Ae ap i + 8Ae equa (and 
consequently the total correlation energy as well) is char- 
acterized by small linear terms and by large quadratic 
terms. The physical explanation is given by a compe- 
tition between the bonds Ti - apical O and Ti - equa- 
torial O. Suppose we distort the crystal in such a way 
that only the internal coordinate r\ increases. The api- 
cal oxygens are a little further apart and the interatomic 
van der Waals-like correlation energy contained within 
Ae ap i decreases in magnitude. At the same time, the po- 
larizing influence of the Ti ions on the charge clouds of 
the apical oxygen anions decreases; this leads to a charge 
displacement from these anions towards their other Ti 
neighbours, with respect to which they are in equatorial 
positions, and, as a result, the electron correlation con- 
tained in |Ae egtta | is getting bigger. Analogous effects (of 
opposite sign) arise with a change of r 2 , with the excep- 
tion that the two equatorial Ti neighbours of an oxygen 
anion are moved simultaneously, which leads to a smaller 
prefactor of Ar 2 , in the expression for Ae equa , as com- 
pared to that of An for Ae ap i- The two discussed types 
of correlation-energy change are opposite and in linear 
approximation almost cancel. 

In conclusion we have shown that electron correlations 
do not change the lattice geometry in a major way but 
are important for other ground-state properties — the 
cohesive energy and bulk compressibility which we ob- 
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tained are close to the experimental values. In view of 
these results, we conclude that the application of the in- 
cremental scheme to rutile yields new insights (in spite 
of the fact that we were forced to go to approximations 
which arc problematic from the methodological point of 
view). In addition we gained valuable experience about 
how to apply this scheme to mixed covalent-ionic systems 
with fractional ion charges. 
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TABLE I. Ab initio results for rutile 



SCF 



SCF b 



SCF present work 



SCF + CCSD(T) 



experiment 



0.0 

co c 

Xo 

B d 



4.559 
3.027 
0.3048 
2.81 



4.555 
3.024 
0.3061 
2.79 



4.529 
3.088 
0.3052 
3.04 



4.548 
2.993 
0.3046 
2.36 f 



4.592 
2.958 
0.3048 
2.39 



Biso c 
AE B C 



3.09 
1.025 



3.08 
1.175 



3.15 
0.834 



2.45 
1.422 



1.470 



Ref. 



Ref. 



data (all-electron calculations) 
data (pseudopotentials on Ti and O) 
" units are A 
d units are Mbar 

c AEb is cohesive energy pes unit cell; in Hartree units, including zero point vibrations with a Debye approximation and a 
Debye temperature of 530 the energies of the free atoms are averaged over J with the appropriate experimental spin-orbit 
contributionsEJ) 
The estimated error bar for B is ±0.06 
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TABLE II. Selected correlation-energy increments in rutile, from CCSD calculations (SCF equilibrium geometry from Ref. 
4, 2nd column in Table I; atomic units). For the free atoms, we obtain correlation energies of -0.1714 H for O, and -0.4040 H 
for Ti, at the CCSD level. 
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a Ti - nearest neighbour apical oxygen (distance - 3.72 a.u.) 

b Ti - nearest neighbour equatorial oxygen (distance - 3.70 a.u.) 

c nearest neighbour oxygen pair (distance - 4.72 a.u.) 

d equatorial oxygen - Ti - apical oxygen (distances: 3.72 and 3.70 a.u.) 

e equatorial oxygen - Ti - equatorial oxygen (distances: 3.70 and 3.70 a.u.; the O-Ti-O angle is 79°.) 

f Ti - O - Ti (distances: 3.72 and 3.70 a.u.) 

g nearest neighbour Ti-Ti pair (distance - 5.71 a.u.) 

h next nearest neighbour oxygen pair (distance - 5.25 a.u.) 



FIG. 1. The rutile structure with oxygen (open circles) and titanium atoms (filled circles). 
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